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Complex kinetic problems are generally 
modeled employing numcrieal integra- 
tion routines. Our kinetics modeling 
program, Aeuchem, has been modified 
to fit rate constants and absorption eo- 
efficients generieally to real or synthe- 
sized "laboratory data" via a least- 
squares iterative procedure written for 
personal computers. To test the model 
and method of analysis the .self- and 
cross-eombination reactions of HO2 and 
CH3O2 radieals of importance in atmo- 
spheric chemistry are examined. These 
radieals as well as other species absorb 
ultraviolet radiation. The resultant ab- 
sorption signal is measured in the labo- 



ratory and compared with a modeled 
signal to obtain the best-fit to various 
kinetic parameters. The modified pro- 
gram generates synthetic data with 
added random noise. An analysis of the 
synthetic data leads to an optimization 
of the experimental design and best- 
values for certain rate constants and 
absorption coefficients. 
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1. Introduction 

Many rate constant determinations, particularly 
for radical-radical reactions, are obtained in com- 
plex kinetic systems in which more than one rate 
process and more than two species concentrations 
are involved. The kinetic processes can be first or 
second order or can be mixed. In such cases, an an- 
alytic relationship involving species concentrations 
and rate constants is often impossible to attain and 
numerical integration is required. If the measure- 
ments involve optical absorption, for example, an 
added complication of contributions to the signal 
from more than one species may arise. There are a 
number of modeling programs available to handle 
such complex systems but most are written for 
main-frame computers and few are capable of 
optimizing the fit of a given model to laboratory 
data by simultaneously adjusting a number of 
parameters [1]. To obtain a "best-fit" of the data to 



the kinetic parameters the least squares procedure 
can be applied. While there is in principle only one 
least squares method [2,3], its application requires 
consideration of the conditional constraints 
relevant to each problem. These are sometimes 
referred to as the condition or adjustment equa- 
tion(s). 

We describe in this paper the application of a 
computer program to a least squares fit that is 
applicable to any mechanism. However, the condi- 
tion equations may be incomplete and the pro- 
cedure for analyzing the fit of the parameters to the 
data should be examined carefully by the user. The 
kinetic modeling computer program, Aeuchem [4] 
is now widely used. It has been modified in the 
present work to: 1) perform "Monte Carlo" simula- 
tions, and 2) to analyze synthetic data (with noise) 
via the least squares method and to yield a best fit 
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to various kinetic parameters within the constraints 
of the experiment. The program can also serve for 
the analysis of laboratory data as well an aid in the 
optimization of experimental design. 

The full program capabilities can be demon- 
strated through the application to the analysis of 
an important laboratory kinetics problem— which 
will be done in this paper. Here, the application of 
the program will be emphasized rather than its de- 
tailed construction. The details of the program will 
be published in near future. However, an inter- 
ested reader wishing to implement methodology 
described here could request the program from 
authors. 

2. Procedure 

The Acuchem modeling program requires an 
input file containing a kinetic mechanism with ini- 
tialized rate constants and initial species concen- 
trations. An output file is generated which lists the 
concentrations of all of the species on a pre- 
selected time grid, and can be displayed as a user- 
directed graph. We have written an additional 
program that converts the Acuchem output file 
into a synthetic "laboratory" data file. The name of 
each species is printed separately and the user can 
enter the contribution of each species in the com- 
posite "experimental" absorption signal. In normal 
practice, the synthetic data would be used to deter- 
mine the sensitivity of the parameter(s) whose val- 
ues are being sought to other measured variables. 

A new analysis program called Acufit has been 
devised. It reads the Acuchem input file, assumed 
to contain the correct mechanism or "model," and 
then either the synthetic "laboratory" or real labo- 
ratory data. This is fitted to the "model" in a pre- 
scribed way. Both rate constants and absorption 
coefficients (representing the contribution of each 
species to the signal) can be adjusted. An iterative 
adjustment procedure is employed. These prob- 
lems are inherently nonlinear and one adjustment 
around initial "first guess values" (by taking only 
first derivatives) is usually not sufficient. Of course 
first guess values are necessary and, as long as they 
are not too divergent from the "correct" values, the 
iterations will converge (if a solution is possible, 
see below). Adjustments are made via solution of 
the "normal equations." Because the variance- 
covariance matrix is readily available, the quality of 
the fit can be assessed. Solving for all possible 
parameters in the system is usually not possible in 
that multiple (degenerate or near-degenerate) so- 
lutions can often fit the data equally well. Programs 



similar to the one described here are not fully auto- 
matic because operational decisions must be made. 
The possible variations in the analysis method are 
explored. 

3. Mechanism (Model) 

The specific kinetic problem we investigate here 
consists of determining the cross-combination rate 
constant for two radicals of atmospheric interest, 
HO2 and CH3O2. Additional details of this system 
can be found elsewhere [5,6]. The experiments are 
complicated by overlapping absorptions from HO2 
and CH3O2, uncertainties in absorption cross sec- 
tions, and the inability to perform the experiment 
under first-order conditions. The objective of this 
paper is not resolving the disagreement concerning 
the correct values for absorption cross sections or 
rate constant for HO2 + CH3O2 but to demonstrate 
the application, capabilities, and limitations of our 
current kinetic modeling program. To do so, we 
have chosen the initial conditions and results of 
earlier studies of this system in our laboratory. Pri- 
marily, because conditions of other studies are not 
very clearly outlined in the literature [5,6]. 

There are three dominant reactions, two combi- 
nation reactions and the cross-combination reac- 
tion. The initial radical concentrations (or rather 
their ratio) are in the range of previously used ex- 
perimental values [5,6]. The units below in Eqs. (1) 
and (2) are cm^ molecule"' s"' for k^, molecules 
cm"^ for concentrations, nanometers for A, and 
10"^" cm^ (base e) for absorption coefficient, &. 

H02 + H02^H202 + 02, ^1 = 1.86x10"'^ 

CH3O2 + CH3O2 -^ products, ^2 = 3.60x 10"'^ 

HO2 + CH3O2 -> CH3O2H -t- O2, k3 = 2.90 X 10- '^ 

[H02]i) + [CH302})= 10.8 X 10'^ A (215) 

[H02]n+[CH302]o= 11.8X 10'^ A (250) 

(1) 

The parameters listed are used to synthesize a 
good surrogate to the "real" laboratory data, which 
can be used for test purposes. The objective here is 
to assess the optimum conditions for obtaining a 
best fit to the cross-combination rate constant fe. 
The rate constants ki and fc and the e's are best 
obtained in single component systems [7,8,9,10] 
containing only HO2 or CH3O2. 

Any complex experiment of this kind entails con- 
ditions that are characteristic of the experimental 
design. Those that apply to an earlier study from 
this laboratory [5] are given in Eq. (2). 
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eH02 X [H02]o + eCHaOs X [CHjOjIo = 
3.73x10"*, A(215) 

eH02X[H02]o+eCH302X[CH302]o = 

3.39x10-", A (250) (2) 

The conditions listed in Eq. (2) arise because, i) 
the sum of the initial concentrations, [H02]i) and 
[CH302](i, results from the titration of known num- 
ber of precursor CI atoms, and ii) for weak absorp- 
tions the measured absorption signal is linear in 
the sum of the product of the absorption coeffi- 
cient and the concentration (and path length) for 
each respective species. The numbers given under 
Eq. (2) must always be consistent with those listed 
under Eq. (1). Also, for convenience, the program 
requires that the ; =0 concentrations be initialized 
and cannot be adjusted by the program. In the 
present work only dimensionless units are 
employed. This does not influence the analysis of 
errors. 

4. Results of Experiments Performed at 
250 nm 

We have applied the Acufit program to the 
model defined by Eqs. (1) and (2) with a view to- 
wards: i) identifying the number of parameters that 
have to be fixed under various noise/signal (n/s) 
conditions; ii) specifying different parametric val- 
ues (in groups) in order to identify those producing 
the least bias in the analysis of ks; iii) fixing all 
parameters with the exception of kj, the cross- 
combination rate constant, in order to provide an 
independent check on the error analysis obtained 
in a previous experimental study; and iv) identify- 
ing initial concentration ratios of reactants which 
produces the least bias on fe given variations of 
known parameters within their expected uncer- 
tainty. Results are presented in tabular and graphic 
form demonstrating the effects of these changes in 
the method of analysis. 

4.1. Fixed Parameters, at "Perfect" Values 

As stated earlier Acufit requires as input: i) a 
synthesized or a real laboratory data file, and ii) 
the Acuchem (model) input data file. In the exam- 
ple we used here the synthesized file was con- 
structed using the values of ^i, ^2, fc, eH02, 
eCH302, eH202, eCH302H, [H02]o, and [CH3O2],,, 
given under Eq. (1) for A =250 nm. The Acuchem 
(model) input file contains all of the information 
listed under Eq. (1), i.e., the mechanism, the initial 



concentrations (always assumed to be known) and 
the rate constants, ki, ^2, and ks taken as initial first 
estimates. Acufit also requires specification of the 
parameters (k's or e's) to be evaluated and their 
initial assignments. Parameters not to be evaluated 
(fixed) are held at their input values. Finally Acufit 
requests the value for the random (normally dis- 
tributed) noise on the data and whether the noise 
distribution is to be always the same (unseeded) or 
randomly variable (seeded) from run to run. In the 
present study the "noise" or a applied to the data 
is taken to be independent of the signal level, com- 
mon for most absorption experiments. 

Since the initial assignments for all of the 
parameters are similar to those used in the con- 
struction of the synthetic "laboratory" file, Acufit 
might return values very close to these, depending 
upon the entered level of noise. In most cases, 
however, the sums of residuals squared either in- 
crease with successive iterations or oscillate such 
that some of the adjusted rate constants or absorp- 
tion coefficients become negative (the program re- 
turns negative ^'s as zero). This means that this 
problem is not sufficiently constrained so that a 
unique solution can be found. If this occurs one or 
more parameters must be fixed at some prescribed 
input value. If that value is the same as used in the 
construction of the synthetic "laboratory data" file, 
it is called a "perfect" value. If it is "fixed" at a 
different value it is called "imperfect." Only when 
"perfect" values are used for the fixed parameters 
will the analysis program return values for the ad- 
justed parameters that are the same as those used 
in the construction of the synthetic file, i.e., within 
the precision given by the variance-covariance 
matrix. The precision, of course, depends on the 
noise level. 

Table 1 was constructed to illustrate variations in 
the method of analysis. Column labels all 
parameters that affect the outcome of the best-fit 
of the data to the model. Values used in creating 
the "synthetic laboratory" file are listed in column 
1. Subsequent columns list results of the analysis 
with error limits obtained from the variance-covari- 
ance matrix. If a cell contains the symbol F, the 
exact value for that parameter as given in column 1 
is employed. If the value is fixed at a value differ- 
ent from that in column 1 it is so designated (for 
example F =0). All columns involve analysis of the 
same "synthetic" laboratory data file. This file 
contains the same noise distribution (at the same 
or different level as specified in the first row). That 
is, in running the program the random number 
generator that produces the distribution is 
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"not seeded." Figure 1 shows the best-fit to the 
synthesized signal at 0.01 and 0.05 nls levels. It is 
noted that in spite of the different intensities, the 
distribution with regard to the noise is the same. 
The best-fit to the synthesized data, in Fig. 1, was 
obtained by fixing all parameters at their "perfect" 



values, (see column 1, Table 1). Most laboratory 
experiments are performed at nls levels between 
these two extremes. In comparing the synthesized 
signal to a real lab signal it is important that the nls 
levels be properly compared particularly if dimen- 
sionless units are employed. 



Table 1. Different analysis fixing various parameters at tlieir "perfect" values and adjusting those remaining 



Para- Values 


n/s = 


n/s = 


nls = 


n/s = 


n/s = 


n/s = 


„/s = 


n/s = 


n;s = 


n/s = 


n/s = 


n/s = 


n/s = 


meters 




0.01 


0.01 


0.01 


0.01 


0.01 


0.01 


0.01 


0.01 


0.05 


0.01 


0.05 


0.01 


0.05 





1 


2" 


3" 


4» 


5" 


6 


7» 


8» 


9 


10 


11 


12 


13 


14 


£lI02 


60 


57.2 ± 
10 


F 


12,8 s 
84. 


161s 
178. 


F 


F 


F 


54.3 s 
11.3 


38. IS 
56.3 


F 


F 


F 


F 


eCIIjOz 


365 


F 


364 s 


380 S 


330 s 


365 S 


F 


F 


367 s 


372 s 


366S 


369. ± 


F 


F 








3.5 


27 


63. 


3.2 


F 


F 


2.7 


10.9 


2.2 


9.3 


F 


F 


eHjOj 


8.3 


F 


F 


F 


F 


F 


F 


f = 0.0 


F 


F 


F 


F 


F 


F 


iClhOiU 


3.98 


F 


F 


F 


F 


F 


F 


f=0.0 


F 


F 


F 


F 


F 


F 


k, 


1.86 


1.25 ± 
0.05 


1.26 s 
1.03 


0.97 s 
1.5 


F 


F 


1..38S 
0.97 


1.31ft 
0.94 


F 


F 


F 


F 


F 


F 


k2 


0.36 


0.34 s 
0.02 


0J4S 
0.03 


F 


0.31s 
0.09 


0.36 s 
0.01 


035 s 
0.02 


0.34 s 
0.02 


F 


F 


F 


F 


F 


F 


k. 


2.90 


3.02 s 


3.01 s 


3.18S 


2.59 S 


3.05 s 


3.05 s 


3.05 s 


3.08 = 


3.82s 


3.08 s 


3.80s 


3.00 S 


3.52s 






0.17 


0.20 


0.23 


0.78 


0.20 


0.12 


0.12 


0.19 


1.12 


0.20 


1.07 


0.10 


063 


lCihO,l, 


0.88 


F 


F 


F 


F 


F 


F 


F 


F 


F 


F 


F 


F 


F 


[HO,!, 


0.30 


F 


F 


F 


F 


F 


F 


F 


F 


F 


F 


F 


F 


F 



"The combination of parameters adjusted in this column cannot be uniquely determined if n/i = 0.05. 




TIIME (ARBITRAIIY UNITS) 



Fig. 1. Best fit to the data when the n/s ratios are 0.01 lower and 0.05 upper (insets). Parameters employed in synthesis 
and analysis are given in column 1, Table I. 
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If a file is produced from a program in which the 
noise distribution is seeded, it is much the same as 
when a given experiment is re-run. Regardless of 
whether the noise distribution is seeded or not the 
error in fitting a parameter should usually fall 
within the error estimate given by the variance- 
covariance matrix. The error analysis may be 
checked by re-running in the seeded mode. Each 
experiment returns a slightly different value for the 
adjusted parameters and the standard deviation 
may be compared with the variance-covariance 
error estimates. Such checks demonstrate that the 
program performs properly. 

For this mechanism or model, the program will 
not be able to adjust all parameters (Appendix A). 
There may be as many as nine parameters. It is 
unusual in such kinetic mechanisms for the ad- 
justed parameters to be independent of one an- 
other. The variance-covariance matrix shows the 
degree of correlation or anti-correlation between 
them. Species concentrations are excluded from 
adjustment, a reasonable programming decision 
that can be relaxed at some later time. For this 
mechanism the maximum number of parameters 
that can be uniquely adjusted is four (Table 1). The 
number depends on the level of noise on the data. 

To obtain the best fit to fc the question arises as 
to whether those parameters which contain consid- 
erable error should be fixed. Table 1 provides an 
answer to the question for at least one case. The 
analysis listed under column 5 shows a poor value 
for eH02, the error being as large as the value of 
the parameter. The error obtained for h is ± 0.78. 
In column 6 the parameter for eH02 is fixed at its 
"perfect" value, eHO2 = 60. As a result a better 
value of A:3 is obtained with an error of ±0.20. In 
this case, provided that eH02 is known well 
enough, it would clearly be better to fix it than 
search for it. If for example, eH02 were fixed at an 
"imperfect" value of 30 while maintaining all of the 
other entries in column 6, /:3 = 3.2 ±0.22, would be 
obtained, a better value than that shown in column 
5. We will explore in detail the effect of fixing 
parameters at "imperfect" values in the next sec- 
tion but emphasize here the well quoted rule that if 
a parameter is obtained with an error estimate 
comparable to its value it may be better to fix that 
parameter. 

Table 1 shows, as expected, that the parameter 
of interest is obtained with better precision when a 
larger number of parameters are fixed. This advan- 
tage is offset by the fact that fixing parameters at 
"imperfect" values introduces a bias in the results. 
Table 1 shows further that the model is very insen- 



sitive to eH202 and eCHjOzH in that setting these 
absorption coefficients equal to has little influ- 
ence on the adjusted parameters (columns 7 and 
8). Since all parameters cannot be adjusted simul- 
taneously, the best strategy is to fix them at their 
"perfect" values. As shown in Table 1 the errors 
given by the variance-covariance matrix for the var- 
ious parameters is closely linear in the nis level 
(columns 9 and 10, as well as 13 and 14). 

The main observations from the analyses listed 
in Table 1 are, i) a number of parameters must be 
fixed, ii) fixing some combination of parameters is 
superior to other combinations. Parameters that 
show the greatest error estimates should preferably 
be fixed, iii) the amenability of adjustment for cer- 
tain parameters is a function of the noise level and 
initial conditions, and finally, iv) if all fixed 
parameters are entered at their "perfect" values 
the variance-covariance matrix error estimates are 
completely reliable. 

4.2. Fixed Parameters at "Non-Perfect" Values 

We now explore the consequence of not entering 
"perfect" values for the fixed parameters in the 
calculations. After all in an actual experiment, we 
do not usually know the exact value for a fixed 
parameter, only an average value with a plus and 
minus uncertainty. The conclusion is that fixing a 
parameter at a "non-perfect" value simply pro- 
duces a bias on the values of the adjusted parame- 
ters yielded by the analysis. The strategy is to fix 
only those parameters causing the least bias on the 
parameter of interest, k^ in the present example. 
When biases are present, estimates of errors (of 
the adjusted parameters) obtained from the vari- 
ance-covariance matrix are likely to be too small. 
Such biases are best investigated by means of the 
present procedure through the use of synthetic 
data. The biases can then be determined, and only 
when they are added to the variance-covariance er- 
ror estimate is there an appreciation of the real er- 
ror in such a system. 

In developing Table 2 a few observations are 
noted. The fixed parameters were generally chosen 
±10% from their "perfect" values, column 1. This 
produces biases that are nearly symmetric. Since 
the program does not provide for adjustment of 
initial concentrations, when either [HO2}) or 
[CH302]i) is fixed at an "imperfect" value it is neces- 
sary to adjust one or both of the absorption coeffi- 
cients, eHOi and eCH302 so that the implicit con- 
straint, [H02]eH02 + [CH302]eCH302 = constant, 
applies. If eH02 and/or eCH302 are not 
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simultaneously adjusted a large bias occurs 
(columns 6 and 7, Table 2). In the next section the 
consequences of fixing eHOz and/or eCH302 at 
"imperfect" values are explored. In such cases the 
constraints of Eq. (2) must be take into account 
explicitly to properly set [HO2I) ar"^ [CH302]tj. 

A relatively small bias in k^ develops regardless 
of which parameter is fixed at an imperfect value. 
If we were to conclude our error analysis at this 
point and analyze real "laboratory" data at the 1% 
nis level we would expect an absolute error of 
about 0.19, Table 1, due entirely to the random 
noise and independent of any bias. If the initial 
concentrations of HO2 and CH3O2 were uncertain 
to the extent of 10% we would expect a small bias 
of ±0.07 (in h) due to CH3O2 (compare fc in 
columns 6 and 8, Table 2, with 2.90) and a bias of 
±0.23 due to HO2 (compare columns 9 and 10, 
Table 2, with 2.90). Assuming the random errors 
and the biases can be combined through vector ad- 
dition the combined error in fc would be: 
0^(^:3) = (0.19)^ + (0.07)- + (0.23)^ = (0.31)-. Addi- 
tionally if /c2 were fixed and its estimated error were 
also at the 10% level, the term (0.44)^ would have to 
be added. If /ci were fixed, a 10% estimated error in 
this parameter, according to Table 2, would require 
adding the additional term (0.135)^ in order to ob- 
tain the combined error in 0-^(^:3) = (0.55)^ (an error 
of 19% in ki). This assumes that the errors in the 
fixed parameters can be properly estimated and 



they are not correlated. An error in both [HO2}) and 
[CH302]ii with no error in eH02 and eCH302 was 
assumed. The error estimates are associated with 
the product of the absorption coefficient and the 
initial concentration for each species. The errors in 
the initial concentration and the e (for each species) 
are in fact distributed. 

The principal conclusions are, i) that when one 
or more parameters are fixed, the bias on /cj (the 
parameter of interest) that each fixed parameter 
produces must be determined, and ii) If the biases 
are small, as in the present case, errors given by the 
variance-covariance matrix are still valid. They en- 
compass the correct or expected results (column 11, 
Table 2). If results fall outside of the error esti- 
mates obtained from the variance-covariance ma- 
trix the composite error can be obtained by vector 
addition of the random error, due to noise, and the 
bias engendered by fixing a parameter. 

4.3. Fixing eHOj and/or eCHjOz at "Non-Perfect" 
Values 

We continue the analysis of Table 2 by fixing the 
two important absorption coefficients at "imper- 
fect" values. This is equivalent to using a erroneous 
value for one or both of these coefficients in the 
analysis of real laboratory data. Here we are both 
compliant and noncompliant with the constraints of 
Eq. (2). Thus, eH02 and eCH302 are individually 



Table 2. Bias produeed on adjusted parameters by fixing scleeted parameters at "non-perfeet" values 



Para- 


Values 


nls = 


nls = 


nls = 


nIs = 


nls = 


nls = 


nls = 


nts = 


nls = 


nls = 


meters 




0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.01 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


eH02 


60 


91.7 


29.4 


49. 


163. 


F 


F 


F 


55. 


66. 


63. ± 
11.0 


eCH302 


365 


354. 


375. 


402. 


330. 


321. 


F 


407. 


F 


F 


F 


eH202 


8.3 


F 


F 


F 


F 


F 


F 


F 


F 


F 


F 


eCH302H 


3.98 


F 


F 


F 


F 


F 


F 


F 


F 


F 


F 


k, 


1.86 


F = 2.Q5 


F = l.dl 


1.18 


2.43 


1.27 


5.44 


2.30 


2.06 


1.52 


0.81 ± 


ki 


0.36 


0.35 


0.37 


F=0.40 


F=0.32 


0.32 


0.44 


0.39 


0.35 


0.37 


0.35 ± 
0.02 


k. 


2.90 


2.76 


3.03 


3.31 


2.43 


2.82 


7.40 


2.96 


2.69 


3.15 


3.28± 
0.19 


[CH3O2I, 


0.88 


F 


F 


F 


F 


F = 1.0 


F = 1.0 


F = 0.79 


F 


F 


F 


[HO2],, 


0.30 


F 


F 


F 


F 


F 


F 


F 


F = 0.33 


F=Q21 


F=Q.ll 
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or simultaneously varied and the initial concentra- 
tions [HO2}) and [CHsOzJo are fixed at values dic- 
tated by Eq. (2), or are fixed at values given in 
column 1 of Table 3. The results show that non-ad- 
herence to the conditional constraints of Eq. (2), 
will not adversely affect the bias on fe (compare 
columns 2 and 3, Table 3) if eHOz is fixed (10%) 
off its "perfect" value. A large bias does occur, 
however, if eCHsOz is fixed ±10% off from its 
"perfect" value and the initial concentrations 
[HOzJo and [CHaOzJo are not initialized in accord 
with Eq. (2) (see results of columns 4 and 5). 
Columns 9 and 10 show a significant bias on fe if 
eHOj is set ±10% and eCHjOj is set ±10% off 
from their perfect values with both offset in the 
same direction. 



4.4 Optimum [HO>]o/[CH302]» 
Obtaining kj at 250 nm 



Condition for 



We explore here how the absolute error in kj 
varies with the ratio [HOzJd/LCHsO:}]. It may be 
guessed that ki is best-fit under conditions where 
the cross-combination product builds up to its 
largest possible value, but at the same time it is 
desired to maximize the absorption signal. The 
noise level in the absorption signal is a function of 
the intensity of the analysis light beam (the signal 
to noise ratio improves linearly with the absorption 



signal). It is not known a priori which of these two 
optimization conditions most influences the mini- 
mization of error in fo or whether other factors also 
contribute. 

Figure 2 shows the result of the error analysis 
obtained by varying the ratio [HOzJ^CHaOzJo. The 
absolute error in ki is obtained by first preparing a 
"synthetic" laboratory data file appropriate to the 
initial concentrations at the specified noise level, 
0--3.66 corresponding to a «/5=0.01, with all of 
the other parameters given under Eq. (1). Then k 
is adjusted (best-fit) fixing all other parameters at 
their "perfect" values. 

The total absorption signal (5), the final cross- 
combination product (Pi==c) and the error in deter- 
mining k3 (E) are plotted in Fig. 2. Since (P) 
maximizes at about a 50/50 composition ratio of 
[CHsOzVCHOz}) and (5) maximizes well to the 
right (side of the figure) a minimum in the error in 
ki would be expected at a composition ratio greater 
than 50/50. In fact a minimum is obtained at about 
40/60. The prediction indicates that factors other 
than the maximum signal and cross-combination 
product enter into the optimization analysis. The 
minimum error in fe occurs over a fairly extensive 
range of initial concentrations. Clearly the opti- 
mization of the composition condition must be 
made by minimizing the error in k^. 



Table 3. Bias produced on adjusted parameters by fixing absorption eoeffieients and eoneentrations at "non-perfeet" values 
while maintaining eonstraints of Eq. (2) 



Para- Values 


n/s = 


,1/S = 


n/s = 


n/s = 


n/s- 


= 


n/s = 


n/s = 


n/s = 


n/s = 


meters 




0.0 


0.0 


0.0 


0.0 


0.0 




0.0 


0.0 


0.0 


0.0 





1 


2 


3 


4 


5 


6 




7 


8 


9 


10 


cHOj 


60 


F = 66. 


F = 66. 


F 


F 


F 




F 


F = 54. 


F=54. 


F = 66. 


eCHjOi 


365 


F 


F 


F = 400. 


F = 400. 


F = 


474. 


F = 328. 


F 


F = 328. 


f = 400. 


cHoOj 


8.3 


F 


F 


F 


F 


F 




F 


F 


F 


F 


eCHsO.H 


3.98 


F 


F 


F 


F 


F 




F 


F 


F 


F 


k, 


1.86 


F 


F 


F 


F 


F 




F 


F 


F 


F 


k2 


0.36 


F 


F 


F 


F 


F 




F 


F 


F 


F 


k3 


2.90 


2.71 


2.96 


2.20 


5.70 


1.8C 


1 


4.90 


3.20 


5.00 


2.00 


[CH3O2],, 


0.88 


F = 0.87 


F 


f =0.79 


F 


F = 


0.65 


F = 1.00 


F = 0.89 


F = 1.01 


F = 0.78 


[HO2J1 


0.30 


F = 0.31 


F 


F = 0.39 


F 


F = 


0.53 


F = 0.18 


F = 0.29 


F = 0.17 


F = 0.40 
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Fig. 2. Absolute error in k-i, E; Absorption signal at 250 nm, 
S; and eross-combination prodiiet plotted vs mole fractions of 
[CHjO^]!, and [HO;}]. Parameters same as in eolumn 1, Table 1; 
and[CH302]ii + [H02]i.= 1.2. 



5. Optimization of [H02]o/[CH302]o to 
Obtain Several Parameters at 215 nm. 

We extend the above analysis by optimizing 
[HO2MCH3O2}) with respect to several of the para- 
meters, taken one at a time. The experiment is at 



215 nm and the appropriate parameters of Eq. (1) 
are used. The results of the analysis are shown in 
Fig. 3. It is noted that the analysis displayed in Fig. 
2 (experiment at A =250 nm) and 3 (A =215 nm) 
can be directly compared with respect to the magni- 
tude of the error on h, common to both figures. 

The optimization conditions for obtaining the 
parameters k] and fo (also eH02 and eCHsOz) are 
obvious but it is instructive to examine the detailed 
shapes of the curves. Clearly for the self reactions 
the error in the determination of the k's is smallest 
when the associated radical concentration is domi- 
nant. The optimum ratio [CH302]o/[H02]() = 50/50 
producing the least error in fe is somewhat different 
from that found in the above case ([CH302]i)/ 
[H02]o = 40/60) probably because at 215 nm €H02 is 
comparable in magnitude to eCH302. Again there is 
no way of guessing the outcome. 

If the ratio of the two initial concentrations were 
to be optimized to obtain the least error in several 
parameters combined, a 50/50 mixture of [H02](( 
and [CHsOzJii would probably be the best compro- 
mise condition. We might expect to be able to solve 
for more than four parameters uniquely, as found 
previously using a nonoptimum [H02]()/[CH302]{) 
ratio. However, it is found that only four parame- 
ters can be assessed, regardless of how small the 
noise level is set. This means that some combination 
of parameters yield fits that are about as good as 
other combinations. The system must be further 
constrained by fixing some of them. Setting [HO2})/ 
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Fig. 3 Absolute errors in fcj, fci, and ky vs mole fraetions of [CHsO;],, and [HO^}, at 215 nm. Rate parameters same as in eolumn 1, 
Table 1; absorption coeffieients (for A =215.), Conditions as in Eq. (1); [CH302]ii + [H02]ii= 1.2. B: Absolute errors in eCH302 and 
eHO: are similarly plotted. 
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[CH302]o significantly off center ([HO2]/ 
[CH302]o = 50/50, in Fig. 3) does have the expected 
effect of further restricting the number of parame- 
ters that can be varied simultaneously. 

The combined error estimate in fe for a 50/50 
mix of [H02]o/[CH302]o (the optimum mixture) is 
recalculated. This is done as previously by deter- 
mining the bias produced in fo by assuming a 10% 
error in both eH02 and eCH302 instead of in 
[H02]() and [CH3O2}). We also include the bias pro- 
duced by 10% errors in k\ and kj, as well as the 
random error produced by a nls=QXil. A some- 
what lower error estimate is obtained, £r=±0.30 
(an error of ±10%) as compared to cr = 0.55 (an 
error of 19%) determined in Sec. 4.2 for nonopti- 
mum mixtures. This error was evaluated by fixing 
all parameters and solving for A3, thereby allowing 
a direct comparison with the estimated error of the 
previous Ref. [5]. Reevaluating the error in A3 ob- 
tained in Ref. [5] where the percent error in both 
eH02 and eCH302 were chosen to be 30%, we find 
that at a 10% level and including the additional 
bias due to 10% errors in ky and kz, we arrive at an 
error of about 20% in A3, similar to the range of 
values obtained here. Even at the 2o- level, this er- 
ror does not completely encompass the range of 
values reported in the literature [9]. Nevertheless, 
a proper accounting for each study would have to 
be performed separately under their unique set of 
boundary conditions, to ascertain whether the re- 
maining differences are attributable to larger er- 
rors in the input parameters or other factors (such 
as secondary chemistry). 

6. Conclusions 

We have shown that the errors returned by the 
variance covariance matrix are reliable when cer- 
tain parameters are adjusted and "perfect" values 
are used for the fixed parameters. Systematic er- 
rors introduced by employing "imperfect" values 
for the fixed parameters can be evaluated (and cor- 
rected for) if a small bias results in the adjusted 
parameter-of-interest (here k{). The use of a 
generic least squares analysis program applicable 
to ail kinetic mechanisms has been described, 
which can be used to optimize the conditions of an 
experiment. Monte carlo simulations are employed 
and the assumption has been made that they are 
reliable surrogates to real "laboratory" data. 

7. Appendix A. Unique Solution 

It is often difficult to ascertain if the parameters 
that are searched for (adjusted) are uniquely 



found. We define a unique fit purely operationally. 
If one or more parameters are initialized (guessed 
at) at different values we should arrive at the same 
ultimate best fit of the parameters (after several 
iterations). If this happens we define the fit to be 
unique. A second criterion involves employing two 
entirely different versions of Acufit. One program 
employs the general method outlined by Went- 
worth [11] which elaborates on Deming's text [2] 
and involves using the SVD method of matrix solu- 
tion [12]. The second program employs the more 
recent Levenberg-Marquardt method [12] and in- 
volves using the Gauss-Jordan matrix solution. 
Both programs yield exactly the same values for the 
adjusted parameters (and variance-covariance ma- 
trix) when a unique solution is possible. They yield 
somewhat different results if a solution is not ob- 
tained; then both result in unrealistic variance-co- 
variance matrices, the first program gives some (or 
all) matrix entries that are far too small, the other 
gives some entries that are much too high. Only 
when both programs yield the same answers do we 
operationally define that "the" solution exists. So 
far we have not encountered a problem where one 
program succeeds in finding a "unique" solution 
while the other does not. 
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